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We report the first direct search for the Kaluza-Klein (KK) modes of Randall- Sundrum gravitons 
using dielectron, dimuon, and diphoton events observed with the D0 detector operating at the 
Fermilab Tevatron pp Collider at y/s = 1.96 TeV. No evidence for resonant production of gravitons 
has been found in the data corresponding to an integrated luminosity of w 260 pb _1 . Lower limits on 
the mass of the first KK mode at the 95% C.L. have been set between 250 and 785 GeV, depending 
on its coupling to SM particles. 

PACS numbers: 13.85.Rm, 04.50.+h, 12.60.-i 



Phcnomcnological models inspired by string theory in 
which there exist additional spatial dimensions have re- 
cently been proposed to remedy some of the defects in the 
standard model (SM). These models may solve the hier- 
archy problem, allow for low-energy gauge coupling unifi- 



cation, and address the issues of flavor and CP- violation. 

The Randall-Sundrum (RS) model Q of extra dimen- 
sions (ED) offers a rigorous solution to a pressing prob- 
lem of the SM — an apparent large hierarchy between 
the Planck scale at which gravity is expected to become 
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strong (Mpi ~ 10 16 TeV) and the electroweak symme- 
try breaking scale (Mew ~ 1 TeV). This is achieved 
through the geometry of a slice of the 5-dimensional Anti- 
deSitter space-time (AdS^), with a curved metric ds 2 = 
cxp(—2kR\ip\)r] llv dx IJ, dx u ~ R 2 dip 2 , where < \<p\ < n is 
the coordinate along the single ED of radius R, k is the 
curvature of the AdS$ space (the warp factor), x^ are 
the convential (3+l)-space-time coordinates, and if* is 
the metric tensor of the Minkowski space-time. A "hid- 
den" (3+l)-dimensional brane (Planck brane) is placed 
at ip = and the second brane (SM brane) is located at 
ip = 7T. Gravity originates on the Planck brane and the 
graviton wave function is exponentially suppressed away 
from the brane along the ED due to the warp factor. 
Consequently, the 0(Mp\) operators on the Planck brane 
yield low-energy effects on the SM brane with a typical 
scale of A,,- = Mpicxp(— knR), where M-p\ = Mpi/^/Sn 
is the reduced Planck mass. Thus, the hierarchy prob- 
lem is solved if A^ ~ 1 TeV, which can be achieved with 
little fine tuning by requiring kR w 10. This is a natural 
solution, as the only fundamental scale in this model is 
Mpi and k ~ R~ x ~ M Pi . 

In the simplest RS model 0,0, the only particles prop- 
agating in the ED are gravitons. Consequently, they ap- 
pear as a Kaluza-Klein (KK) tower of massive excitations 
from the point of view of the SM brane and can be reso- 
nantly produced in pp collisions. The masses and widths 
of the KK-excitations are related to the parameters of the 
RS model. The zeroth KK mode (G^) remains massless 
and couples to the SM fields with gravitational strength, 
1/Mpi, while the excited modes couple with the strength 
of 1/A W . The excited modes can decay into fermion- 
antifermion or diboson pairs, leading to the characteris- 
tic resonance structure in their invariant mass spectrum. 
In this Letter we report on a search for the first excited 
mode of the KK graviton, G^ 1 ', in the dielectron, dimuon, 
and diphoton final states. Since the graviton has spin 2, 
its decay products arc cither found in the s-wave (dipho- 
tons) or p-wave (dileptons). This leads to the branching 
fraction of the graviton decay in a single dilepton channel 
(U) to be half that of diphotons. 

Phenomenologically, it is convenient to express the two 
RS parameters k and R in terms of two direct observ- 
ables: the mass of the first excited mode of the graviton, 
Mi, and the dimensionless coupling to the SM fields, 
k/Mpi, which governs both graviton production cross 
section (~ (fc/Mpi) 2 ) and the width of the graviton res- 
onance. The theoretically preferred range for M\ is be- 
tween a few hundred GeV and a few TeV, while fc/Mpi 
is expected to be between 0.01 and 0.1. Larger values of 
the coupling would render the theory non-pcrturbative, 
while smaller would require an undesirably large amount 
of fine-tuning. Indirect limits on RS model parameters 
come from precision electroweak data (dominated by the 
S parameter) . There have been no dedicated searches 
for RS gravitons to date. 



We used the D0 detector operating at the Fermilab 
Tevatron pp Collider at y/s — 1.96 TeV with approxi- 
mately 246 pb _1 of data accumulated with dimuon trig- 
gers and 275 pb -1 of data collected with single or di- 
electromagnetic (EM) triggers for this search. To max- 
imize the reconstruction efficiency for dielectrons and 
diphotons, we did not use tracking confirmation and com- 
bined these two channels in a single, calorimeter-based, 
"diEM" channel. The detector, data acquisition system, 
and triggering are detailed elsewhere Q. 

Offline, we required EM objects to have transverse en- 
ergy Et > 25 GeV, be isolated in the calorimeter and 
tracker, have significant fraction of their energy deposited 
in the EM calorimeter, and have their EM shower shape 
consistent with that expected for an electron. We ac- 
cepted EM objects in the central (\r]d\ < 1.1) H and 
forward (1.5 < \rjd\ < 2.4) regions of the calorimeter, but 
required at least one of them to be central. 

The overall efficiency per electron was determined us- 
ing Z — > ee events, and is (91 ± 2)% in the central and 
(82±2)% in the forward regions. The efficiency is uniform 
in Et and r/d, with the exception of the region close to 
the boundaries between the central and forward calorime- 
ters, 1.0 < \r]d\ < 1-1, where it drops by a factor of two. 
Monte Carlo (MC) simulations show that the efficiency 
per photon is 5% lower than that per electron. An addi- 
tional inefficiency of 7% per event arises from the trigger, 
EM objects lost in azimuthal cracks between the central 
calorimeter modules or overlaps with jets in the events. 

Muons were identified in the muon spectrometer (cov- 
ering \r)d\ 2.0) and were required to have a matching 
track in the central tracking detector, transverse momen- 
tum pt > 15 GeV, be isolated, and pass additional hit 
and track quality requirements. Since the muon momen- 
tum resolution degrades rapidly at high px, high-mass 
dimuon events sometimes have the momentum of one of 
the muons misreconstructed. To remedy this and reduce 
non-Gaussian tails in the invariant mass resolution, we 
assigned both muons the same value of transverse mo- 
mentum, based on the weighted average (in 1/pt) of their 
individual pr's. This results in w 30% decrease in the 
RMS of the invariant mass distribution at the cost of a 
modest (« 1%) decrease in the invariant mass resolution. 
To reduce cosmic ray background, muon arrival times in 
the muon detector were required to be consistent with 
that for particles originating from beam collisions. The 
two muons in the event were not required to have oppo- 
site signs, as the sign determination efficiency degrades 
fast at high px- The overall selection efficiency per muon 
is (80 ± 4)%, as determined using Z — > events. 

The above requirements result in 22,786 (17,128) diEM 
(dimuon) events used in the analysis. The main back- 
ground to the RS graviton signal is Drell-Yan (DY) pro- 
duction in the dielectron and dimuon decay channels 
and direct diphoton production in the diphoton channel. 
These backgrounds were estimated using the leading or- 
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FIG. 1: Invariant mass spectrum in the a) diEM and b) dimuon channels. The points with error bars are data and the solid 
line is the overall background (dark shading in a) represents instrumental background). Also shown with an open histogram is 
the signal from an RS graviton with mass Mi — 300 GeV and coupling fc/Mpi = 0.05. 



der (LO) MC generator of Ref. |5j, augmented with a 
parametric simulation of the D0 detector [|| . The sim- 
ulation accounts for the calorimeter and tracker resolu- 
tions, primary vertex position, detector acceptance, and 
pr of the It or 77 system. The CTEQ5L set of parton 
distribution functions (PDF) was used in the simulations. 

In the diEM channel, an additional instrumental back- 
ground arises from QCD multijet and direct photon 
events, with one or more jets reconstructed as EM ob- 
jects. This background was estimated from the data, 
by inverting the shower shape quality requirement, with 
the absolute normalization obtained in the Z-boson mass 
peak. The only other background is due to tt production 
with r's decaying via the electron or muon channel and 
is negligible at invariant masses above the Z-peak where 
the search is performed. Figure ^ shows the invariant 
mass spectrum in the diEM and dimuon channels and 
demostrates good agreement between the data and ex- 
pected background. 

The RS graviton signal was simulated with the 
pythia MC event generator with the CTEQ5L PDF, 
followed by the parametric simulation of the D0 detec- 
tor. The LO PYTHIA cross section was scaled by a con- 
stant if -factor of 1.34 to account for next-to-LO (NLO) 
effects, recently calculated 9] for graviton exchange and 
shown to be similar to those for SM DY production. We 
set limits on the ratio of the graviton production cross 
section and the next-to-NLO (NNLO) pp — ► Z — ► ee 
cross section of 254 ± 10 pb [l0(. Since the Z-peak is 
found in the candidate sample, this approach allows for 
in situ calibration and reduces the overall systematic un- 
certainty. We quote the limits on production of gravitons 
in terms of the absolute cross section, which is obtained 
by multiplying the limits on the ratio by 254 pb. 

A simulated signal is shown in Fig. ^ for Mi — 300 
GeV and k/ Mp\ = 0.05. Since the muon momentum was 
measured in the tracker, while the EM energy was deter- 



mined from the calorimeter, the difference in resolutions 
for the two detectors explains that the mass resolution 
in these two channels is so different. We used a conser- 
vative estimate of the muon momentum and EM energy 
smearing parameters in the detector response simulation 
by attributing the measured width of the Z boson to 
the constant resolution term, which dominates at high 
masses. This choice leads to a somewhat broader than 
expected reconstructed signal and to conservative limits 
on signal cross section. 

To set limits on graviton production, we performed 
analyses in a series of overlapping windows correspond- 
ing to different graviton masses. The width and position 
of the windows were optimized to give the highest sig- 
nal sensitivity via a modified method of Ref. |ll| , which 
takes into account Gaussian fluctuations of an exponen- 
tially falling background. For the diEM channel at high 
masses (> 300 GeV), the background is small so a sym- 
metric window with the width set to six times the de- 
tector resolution was used to maximize the sensitivity. 
Since the muon momentum resolution effects on the in- 
variant mass are very asymmetric and result in a long 
high-mass tail (see Fig.nj, only the lower mass bound is 
used in the dimuon channel windows. Since the internal 
graviton width is negligible compared to the instrumen- 
tal resolution in the range of Mi and fc/Mpiwe studied, 
the window size did not depend on k/Mp\. The overall 
geometrical acceptance for the signal in the diEM chan- 
nel varies between 45% and 62%, depending on the mass 
point. In the dimuon channel, the corresponding varia- 
tion is between 55% and 67%. 

The results of the counting experiments are listed in 
Table I. As the number of events in each window is con- 
sistent with the expected background (the significance 
of an upward fluctuation in the diEM channel at 400 
GeV is < 2 standard deviations), we set limits on the 
graviton production cross section. The limits were set 
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TABLE I: Counting experiments and 95% C.L. upper limits (in fb) on cr(pp —* C?' 1 ' — > it). All masses are expressed in GeV. 
Graviton DiEM Channel Dimuon Channel Combined 



Mass Window Background Data Limit Sensitivity Window Background Data Limit Sensitivity Limit 
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independently in the two channels using a Bayesian tech- 
nique |12| with a flat prior for the signal and systematic 
uncertainties on signal and background taken into ac- 
count. In the diEM channel, the signal uncertainty is 
9%, dominated by the mass-dependence of the EM effi- 
ciency (5%), acceptance calculation (5%), and difference 
between the photon and electron efficiencies (5%). In 
the dimuon channel, the signal uncertainty is 8%, dom- 
inated by the acceptance uncertainty (7%). The com- 
mon source of systematics for both channels is the 4% 
Z boson NNLO cross section uncertainty. The SM back- 
ground uncertainty is 9-12% and dominated by the K- 
factor mass dependence (5%), efficiency determination 
(7% in the diEM and 5% in the dimuon channels), mod- 
eling of the momentum smearing (6%, dimuon channel), 
and the PDF dependence (5%). The uncertainty on the 
instrumental background in the diEM channel is domi- 
nated by the low statistics of the background sample at 
high masses. 

The 95% C.L. upper limits on o{pf> G*W -> it) 
are listed in Table I and shown in Fig. [3 Also shown 
is the expected sensitivity of the search in each channel, 
defined as an average limit expected given the Poisson 
distribution of the background around its mean. We fur- 
ther combined the diEM and dimuon limits after taking 
into account common systematic uncertainties. The com- 
bined limits are very close to the diEM limits (and in fact 
are slightly less restrictive due to the overall small excess 
of observed events in the dimuon channel). We translate 
the limits on the cross section times branching fraction 
into limits on the RS model parameters M\ and k/M-p\, 
as shown in Fig. We did not include an uncertainty on 
the signal cross section related to the PDF and higher- 
order QCD effects. Assuming that it is similar to that 
for DY production (« 10%) increases the cross section 




Graviton Mass (GeV) 

FIG. 2: The 95% C.L. upper limits on a(pp -> G* (1) -> it), 
as a function of the graviton mass. The upper (lower) solid 
line with points corresponds to the dimuon (diEM) channel. 
The dashed lines with points represent the expected limits. 
Also shown with a series of smooth lines the production cross 
sections for various values of fc/Mpi between 0.01 and 0.1. 

limits by 2.5%. This translates into a negligible (« 1%) 
fractional change in our limits on fc/Mpi for any graviton 
mass. 

To conclude, we have performed the first dedicated 
search for Randall-Sundrum gravitons in the dielectron, 
dimuon, and diphoton channels using 246-275 pb _1 of 
data collected by the D0 experiment in the Run II of 
the Fermilab Tevatron Collider. We see no evidence for 
resonant production of the first Kaluza-Klein mode of 
the graviton and set the most restrictive limits on the 
RS model parameters to date. Graviton masses up to 
785 (250) GeV are excluded for k/M Pl of 0.1 (0.01). 

We thank the staffs at Fermilab and collaborating in- 
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FIG. 3: 95% C.L. exclusion limits on the RS model parame- 
ters Mi and k/M-p\. The light-shaded area has been excluded 
in the dimuon channel; the medium-shaded area shows the ex- 
tension of the limits obtained in the diEM channel; the dotted 
line corresponds to the combination of the two channels. The 
area below the dashed-dotted line is excluded from the pre- 
cision electroweak data (see Ref. 0). The dark shaded area 
in the lower right-hand corner corresponds to > 10 TeV, 
which requires a significant amount of fine-tuning. 
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